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The mass m,,^ and total width r^°j of the r?c meson have been measured in two-photon interactions 
at the SLAC e'^e~ asymmetric B-Factory with the BABMi detector. With a sample of approximately 
2500 reconstructed 77c K'^K'^-k^ decays in 88 fb~^ of data, the results arc m^^ = 2982.5 ± 1.1 
(stat) ± 0.9 (syst) MeV/c^ and r;'^^^ = 34.3 ± 2.3 (stat) ± 0.9 (syst) MeV/c^ Using the same decay 
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mode, a second resonance with 112 ± 24 events is observed with a mass of 3630.8 ± 3.4 (stat) ± 
1.0 (syst) MeV/c^ and width of 17.0 ± 8.3 (stat) ± 2.5 (syst) MeV/c^. This observation is consistent 
with expectations for the ric{2S) state. 

PACS numbers: 14.40.Gx, 13.25.Gv 



The mass and width of the rjc meson {J^^ = '"), the 
lowest lying state of charmonium, are not as well estab- 
lished as those of the J/ip meson. The world average [3| of 
the total width is T^^^ = IG.Otg j MeV/c^, with individ- 
ual measurements ranging from 7MeV/c^ to 27MeV/c^ 
with large errors. Recent measurements extend from 
17MeV/c2 to 29MeV/c2. 

A radial excitation of the rjc, the r]c{2S) state, is pre- 
dicted by heavy quark potential models to lie below the 
DD threshold 0. The hyperfine separations {ric,J/4') 
and {ric{2S),ip{2S)) are directly related to the spin-spin 
interaction. These calculations predict the mass splitting 
'™i/'(2S) ~ ™??c(2S) to be in the range 42-103 MeV/c^. The 
Crystal Ball Collaboration observed a peak at 91 ± 
5 MeV, in the inclusive photon spectrum oi^pi^S) decays, 
with a width T < 8 MeV (95% confidence level). This 
peak was considered most likely to be due to il^{2S) 
r]c{2S)^, with the 'rjc{2S) state having a mass of 3594 
± 5MeV/c^. The Belle Collaboration recently reported 
signals attributed to the 77c (25') state, but with substan- 
tially higher masses: for the K'^K~Tr~^ mass distribution 
in exclusive B — > KK'^K~Tr~^ decays [3, they measured 
3654 ± 6 (stat) ± 8 (syst) MeV/c^ and T < 55 MeV/c^ 
(90% confidence level); from a signal observed in the in- 
clusive J/ip spectrum in e"*"e~ annihilation 0, they mea- 
sured 3622 ± 12 MeV/ c^. This state was unsuccessfully 
searched for in ^ X ^ 77 j3| and 77 hadrons 
However, an estimate of the two-photon production 
rate of the rjc(2S) suggested that this meson could be 
identified in the current e+e^ i3-factories. 

In this analysis we measure the masses and widths of 
the rjc and of a state interpreted as the ric{2S) meson, 
by reconstructing 77 ^ — > K''^K^tt^ (X!? ^ tt+tt") 
events in the BABAR detector at the PEP-II energy- 
asymmetric e"'"e~ storage ring at SLAC. The data sample 
was collected both on and slightly below the T(4iS') res- 
onance, and corresponds to an integrated luminosity of 
88 fb-i. 

The BABAR detector is described in detail in refer- 
ence The momenta of charged particles are mea- 
sured and their trajectories reconstructed with two detec- 
tor systems located in a 1.5 T solenoidal magnetic field: a 
five-layer, double-sided silicon strip vertex tracker and a 
40-layer drift chamber. Both devices provide dE/dx mea- 
surement. Charged particle identification is provided by 
a detector of internally reflected Cherenkov light, com- 
plemented by the dE/dx measurement. The energies of 
electrons and photons are measured in a calorimeter con- 
sisting of 6580 CsI(Tl) crystals. 

The mesons are formed by the interaction of two vir- 



tual photons. Since the e"*" and e~ scatter through too 
small an angle to be detected, the two photons are quasi- 
real and nearly aligned with the incident beams. A pre- 
selected sample comprises events having four charged 
tracks with a net zero charge and with total laboratory 
energy less than 9 GeV. This removes most events coming 
from B meson decays. 

A further selection of events is aimed at maximizing 
the ratio S/y/{S + B), where S is the signal and B the 
background, both taken within a ±50 MeV/ window 
around the r]c peak. Events with total transverse momen- 
tum in the center-of-mass greater than 1.05 GeV/c or with 
total energy of neutral particles greater than 0.7 GeV are 
rejected. In order to identify rjc K^K^tt^ events, de- 
cays with one K'^ —>■ tt+tt^ candidate that lies within the 
window 0.482 < M{K'^) < 0.512 GeV/c^ are selected. Of 
the two remaining tracks, we require that one and only 
one be identifled as a kaon; the other one is assumed to 
be a pion. The angle between the momentum and its 
flight path, as determined by the i^" and K^ir^ vertices, 
is required to be small {cos9{K'^) > 0.992). Finally the 
K°K^Tr^ vertex is fitted, with the K° mass constrained 
to the world average value Q . 

The resulting K'^K^n^ mass spectrum is shown in 
Fig. n with a large peak at the rjc mass and a smaller 
peak at the J/tp mass. Although the J/ip cannot be pro- 
duced in two-photon fusion, it is expected to be pro- 
duced with hard photon emission by initial state radia- 
tion (ISR). The boost of the asymmetric collider brings 
the decay products of J/ip mesons travelling in the back- 
ward direction into the acceptance of the detector. 




2.5 3 3.5 



M(Kg K 71) (GeV/c^) 

FIG. 1: The {KsK^tv"^) mass spectrum fitted (solid line) to 
r]c + J/tp + background, as explained in the text. The dashed 
line shows the background component of this fit. 

A thorough understanding of the experimental reso- 
lution is essential to determine the width of the r]c me- 
son. The resolution for the J/4> can be inferred from 
data since its natural width is negligible. This is not 
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the case for the rjc, which has a natural width some- 
what larger than the detector resolution. To help de- 
termine the resolution for the r]c, Monte Carlo calcula- 
tions were performed. The generator jl^l used to sim- 
ulate 'y^ ^ rjc ^ K'^K^TT^ events applies the formal- 
ism of Budnev et al. to calculate the cross-section 
for the process e~^e~ —^ e~^e~^j —^ e^e~r]c- Monte 
Carlo calculations were also performed to generate J/ip 
events produced in e+e^ annihilation with initial state 
radiation. Both rjc and J/ip were assumed to decay into 
K'^K^tt^ with a phase-space distribution. In the Monte 
Carlo simulation, the reconstructed r/c and J/^ masses 
are both shifted by —1.1 MeV/c^ (with statistical errors of 
0.1 MeV/c^ and 0.2 MeV/c^ respectively) from their gener- 
ated values. This bias does not affect the mass difference 
m j/^ — m^^ . The mass resolution is estimated by fitting 
the distribution of the difference between reconstructed 
mass and generated mass to a Gaussian function. Its 
standard deviation is found to be 7.3 ± O.lMeV/c^ for 
the r]c and 8.1 ± 0.2MeV/c2 for the J/?/'. 

To determine the mass and width of the rjc, an un- 
binned maximum hkelihood fit to the K'^K^ TT^ mass 
spectrum for masses between 2.5 and 3.5GeV/c^ is per- 
formed. The rjc is represented by a Breit-Wigner func- 
tion (r/2)V((VF - TO^J2 + (r/2)2), with W the invari- 
ant KgK^n^ mass, convolved with a Gaussian resolu- 
tion function. The J/ip peak is fitted with a Gaussian 
function. The background is represented by an exponen- 
tial function of W, Aexp{—XW). The free parameters 
of the fits are the J/tp mass mj/^ , the mass difference 
nT-j/ip — TTT'Tja, the rjc width T^^^, the J/ip resolution aj/^ , 
the coefficients A and A of the background, and the num- 
bers of events in the rjc and J/ip peaks. The resolution 
ar,^ of the r]c peak is constrained to a value 0.8MeV/c^ 
lower than the J/ip resolution, as indicated by the Monte 
Carlo simulation. The results of the fit are: mj/^ = 
3093.6 ± 0.8MeV/c2, mj/^ - m,,^ = 114.4 ± 1.1 MeV/c^, 
CTj/^ = 7.6 ± 0.8MeV/c2, r^'^^^ = 34.3 ± 2.3MeV/c2. The 
numbers of rjc and J/tp events are respectively 2547 ± 90 
and 358 ± 33. 

The mass resolution found for the J/ifj is 0.5 ± 
0.8MeV/c^ lower than the Monte Carlo prediction, but 
consistent with it. To evaluate the systematic uncer- 
tainty affecting the rjc width, the conditions of the fit are 
varied as shown in Table 1 . When aj^^ and cr^^ are fixed 
to the values obtained in the Monte Carlo simulation 
(second row of Table 1), the width of the rjc changes by 
0.6 MeV/c^. We take this value as an estimate of the sys- 
tematic uncertainty associated with the uncertainty on 
the rjc resolution. The value of T^^^ changes by 0.4 MeV/c^ 
on average when the mass range of the fit is varied from 
2.4-3.6 GeV/c^ to 2.7-3.3 GeV/c^. This gives an estimate 
of the systematic uncertainty associated with the choice 
of the mass range of the fit. By varying the event selec- 
tion parameters, we estimate that the systematic uncer- 
tainty associated with the event selection is 0.5MeV/c^. 



The total systematic uncertainty on the rjc width is then 
0.9MeV/c^. The final value of the rjc width is: 

r1;^t = 34.3 ± 2.3{stat) ± 0.9{syst) MeV/c^ 

TABLE I: Results of unbinned maximum likelihood fits to 
the rjc and J/tp mass spectra. The resolutions of the J/tp 
and rjc peaks are respectively ctj/^ and cr,,^. The first row 
presents the nominal fit, and the succeeding rows are used for 
systematic studies of the rjc width. "MC" denotes results of 
Monte Carlo simulations. 



mass range 
MeV/c^ 


p'7c 

tot 

MeV/c^ 


MeV/c^ 


MeV/c^ 


2.5-3.5 
2.5-3.5 


34.3 ± 2.3 
33.7 ± 2.0 


7.6 ± 0.8 
8.1 (MC) 


7.3 (MC) 


2.4-3.6 

2.6- 3.4 

2.7- 3.3 


33.7 ± 2.3 
34.4 ± 2.3 
34.7 ± 2.4 


7.6 ± 0.8 

7.7 ± 0.9 
7.7 ± 0.8 


aj/^-0.8 
cTj/V- -0.8 



The rjc mass is 2982.5 ± 1.1 (stat) MeV/c^, obtained 
by subtracting 114.4 MeV/c^ from the current world aver- 
age value of the J/tp mass JlJ . The rjc and J/tp masses are 
unchanged by the alternative fits listed in Table 1. We 
estimate that the systematic uncertainty on to j/^ — m^^ , 
associated with the event selection, is 0.8MeV/c^. After 
correction for the — l.lMeV/c^ shift seen in simulation, 
as mentioned above, the J/tp mass is still shifted by an 
additional — 2.2Me\Yc^ relative to the well established 
world average value m . Because J/ip events and rjc events 
populate different regions of detector acceptance, as illus- 
trated in Fig. 121 for final-state pions, a shift that applies 
to the J/tp may not entirely apply to the rjc due to pos- 
sible imperfections in the detector modelling. When one 
selects rjc events with decay particles going backward, as 
is the case for the J/tp , the rjc peak shifts by 0.5 MeV/c^, 
which we take as a contribution to the systematic uncer- 
tainty. The final value of the rjc mass is then: 

m,,^ = 2982.5 ± l.l(stat) ± Q.9{syst) MeV/c^. 

The peak at 3.63GeV/c^ in the K^Ktt mass spectrum 
(Fig. m may be the expected rjc{2S) state. In order to 
optimize the significance of the signal, a new event selec- 
tion is performed that maximizes the ratio S/^/B. This 
is appropriate in place of S/ \/S + B because we need to 
establish the significance of the peak without bias from 
assumptions about how much signal to expect, and in 
any case the branching fraction and 77 width needed 
for such a prediction are unknown. For S we take the 
signal as generated from Monte Carlo simulation and 
B is the background estimated from the average of the 
rjci2S) sidebands 3.30-3.48 GeV/c^ and 3.78-3.96 GeV/c^ 
of the data. Events are required to have the total en- 
ergy deposited by neutral particles less than 0.25 GeV 
and cos 6{K'^) > 0.995. The other optimized cuts are the 
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FIG. 2: Angular distributions of pious from the decays of 
J/tp , r]c aud rjc{2S), iu the laboratory frame (Qt^ is the pion 
polar augle). The backgrouuds determined from sidebands 
have been subtracted. 
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FIG. 3: The K°sK'^ mass spectrum with event selection 
optimized for the r]c(2S) as described in the text. The solid 
curve is the fit with the ryc(2S) resonance shape being repre- 
sented by a Breit-Wigner function convolved with a Gaussian 
resolution function. The dashed curve shows the background 
component of this fit. 



same as for the r^c- The resulting mass spectrum is shown 
in Fig. El 

The mass resolution determined from the Monte Carlo 
simulation is 9.2MeV/c^ and the reconstructed mass is 
0.4MeV/c^ lower than the generated mass. Since the res- 
olution for the J/ip was found to be 0.5 ± 0.8MeV/c^ 
lower in the data than in the Monte Carlo simulation, we 
assume that the resolution for ric{2S) is also 0.5MeV/c^ 
lower in the data, with an uncertainty of 0.8MeV/c^. 
The KgK^TT^ mass spectrum is then fitted between 
3.3GeV/c^ and 4.0GeV/c^, the ric{2S) resonance shape 
being represented by a Breit-Wigner function convolved 
with a Gaussian resolution function with standard de- 
viation 8.7MeV/c^. The background is fitted with an 
exponential shape. The fit results in 112 ± 24 events in 
the rjc{2S) peak. The significance of this signal is charac- 
terized by the quantity ^2 x log Cmax I — 4.9, where 
Cmax and £o are respectively the likelihoods for the fits 
with and without the r]c{2S) peak. 



The 



mass difference is found to be 



534.6 ± 3.4 (stat) MeV/c^. Taking into account the shifts 
from generated to reconstructed masses of — l.lMeV/c^ 
for the J/ip and -0.4MeV/c^ for the 7]c{2S), as found in 



the Monte Carlo simulation, this mass difference becomes 
533.9 MeV/c^. The r]c{2S) mass is then m^^(2S) = wj/v- 
+ 533.9 3630.8 ± 3.4 (stat) MeV/c^. The measured 
total width is 17.0 ± 8.3 (stat) MeV/c^ The resolution 
uncertainty of 0.8MeV/c^ results in a systematic uncer- 
tainty of O.lMeV/c^ on the ■qc{2S) mass and 2.0MeV/c2 
on its total width. When the mass range for the fit is var- 
ied to 3.2-4.1 or 3.4-3.9 GeV/c^, the 1)^(28) mass varies 
by 0.2MeV/c^ whereas its width varies by 1.2MeV/c^ on 
average. The 0.5 MeV/c^ uncertainty on the —2.2 MeV/c^ 
shift observed for the measured J/tp mass relative to the 
world average value is taken as a systematic uncertainty 
on the r]c{2S) mass. Based on the upper limit for the 
branching fraction ip{2S) K~^K~tt^ [iLwe estimate 
that tpi2S) (with a mass of 3.686 GeV/cn3|) could con- 
tribute up to 5 events to the spectrum of Fig. 3. Allow- 
ing for this reduces the rjc{2S) width by 0.7 MeV, which 
we take as a systematic uncertainty, whereas the r]c{2S) 
mass varies by about 0.1 MeV/c^. The systematic uncer- 
tainties associated with the event selection are taken to 
be the same as for the rjc, 0.8 MeV/c^ for the r]c{2S) mass 
and 0.5MeV/c^ for its total width. Adding all systematic 
uncertainties in quadrature, the final results are: 

TOr,<,(2S) = 3630.8 ± 3A{stat) ± 1.0{syst) MeV/c^ 



= 17.0 ± 8.3{stat) ± 2.5{syst) MeV/c^ 




0.2 0.4 0.6 0.8 

P.rtot(c.m.) (GeV/c) 

FIG. 4: Total transverse momentum in the center-of-mass. 
The hatched solid line is the result of the two-photon Monte 
Carlo simulation for the rjc{2S) state, normalized to the 
data. The data are events in the 3.60-3.66 GeV/c'^ mass re- 
gion; the background determined from mass sidebands 3.30- 
3.48 GeV/c^ and 3.78-3.96 GeV/c^ has been subtracted. 

While we have not measured the quantum numbers 
of the state at 3630.8 MeV/c^, demonstrating that it is 
formed from the fusion of two quasi-real photons would 
at least restrict the possibilities. Such a process can occur 
only if C=+, and = 0^ (0"*" is excluded by the final 
state), 2=*=, 3+, 4^, .... Other combinations would be pos- 
sible if production were via an ISR process, or if at least 
one of the two photons in two-photon fusion were highly 
virtual. However ISR is excluded as the source, because 
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the decay products of this state have angular distribu- 
tions concentrated in the forward hemisphere, hke the rjc^ 
in contrast to the J/ijj for which the decay products peak 
in the backward direction. This is illustrated in Fig. |21 
Moreover the distribution of the total transverse momen- 
tum (Fig. ^ is peaked at 0, characteristic of quasi-real 
photons, and this excludes spin-one production. Thus 
the evidence supports the state having quantum num- 
bers J^^ ^ 0"+ or J > 2. But J > 2 is disfavored for a 
charmonium state of such low mass, which suggests that 
the state has the quantum numbers of the 77c (25'). 

In summary, we have measured the mass difference be- 
tween the J/tp and the rjc and the total width of the ?7c, 
using 2547 ± 90 events of 77 ^ ?7c ^ X°X±7r=F and 358 
± 33 J/i}) -> K^K^TT^ events, selected with the BABAR 
detector. 

A state which could be the expected r]c{2S) was also 
observed in the K^K^tt^ decay mode, with 112 ± 24 
events, and its mass and total width measured. The 
measured mass is significantly different from the mass of 
the state reported by the Crystal Ball Collaboration 0, 
but consistent with the measurements of the Belle Col- 
laboration 1^ . We have presented evidence that this 
state is produced via the fusion of two quasi-real photons, 
which suggests that its quantum numbers are those of the 
r]c{2S). The deduced mass splitting m^(2S) — Wr;<,(2S) = 
55.2 ± 4.0MeV/c^ is consistent with theoretical expecta- 
tions. 
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